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ABSTRACT 

Bright X-ray flares are routinely detected by the Swift satellite during the early afterglow of 
gamma-ray bursts, when the explosion ejecta drives a blast wave into the external medium. We 
suggest that the flares are produced as the reverse shock propagates into the tail of the ejecta. 
The ejecta is expected to contain a few dense shells formed at an earlier stage of the explosion. 
We show an example of how such dense shells form and describe how the reverse shock 
interacts with them. A new reflected shock is generated in this interaction, which produces 
a short-lived X-ray flare. The model provides a natural explanation for the main observed 
features of the X-ray flares — the fast rise, the steep power-law decline, and the characteristic 
peak duration Af/f = (0.1 - 0.3). 
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1 INTRODUCTION 

Abo ut 30% of GRBs show X-ray flares during the i r early afterglow 
(e.fi. iBurrows et all l2005bl : Ipalcone et alj 1 20071 : IChincarini et al.l 
bold) . Sometimes they are acc ompanied by significant flux in¬ 
crease in the optical band (e.g. lUi et alj|2012h . The X-ray flares 
are characterized by a fast rise of luminosity (typically by a 
factor ~ 10) to a sha rp peak followed by a power-law decay 
dChincarini et al.l[20071) . The ratio of the characteristic temporal 
width of the peak Atf to its time of occurrence since the beginning 
of the GRB is typically Atf/tf ~ 0.1 - 0.3. At the end of the flare, 
the X-ray flux resumes the underlying smooth decay of the early 
afterglow, with no apparent flux increment left. 

Besides the X-ray flares, the early afterglow shows other 
puzzling features, such as plateaus and sudden steep decrease 
in the observed luminosity. These fea tures are not explained 
by the standard forward shock model jMeszaros & ReesI Il997l : 
ISari et al]|l998h . and it was proposed that the observed afterglow 
is produced by a long-lived reverse shock ins ide the GRB ejecta 
dUhm & Beloborodo\l2007l : lGenet et al.11200'^ . The dynamics and 
emission of the reverse shock are particularly sensitive to the struc¬ 
ture of the ejecta (the distribution of its Lorentz factor, density, 
and magnetic fields) which may explain the rich phenomenology 
of the early afterglow jHascoet et al.ll201 l1 . l2012l : luhm et al.ll2012 : 
iHascoet et al.ll20149) . 

The fast rise and relatively short duration of the X-ray flares 
motivated several authors to invoke late activity of the central en¬ 
gine and associate the flares with the jet produced by this activ¬ 
ity — a scaled-down version of the prompt GRB emission (e.g. 
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iBurrows et al. I l2005al: iFan & Weil IIOO^: IZhang et alllloo^ . This 
possibility requires a mechanism that keeps the central engine ac¬ 
tive for about one day, as some flares are observed as late as a day 
after the GRB explosion. It also requires the engine to be variable 
in a special way on a timescale Aff ~ (0.1 - 0.3)ff, in contrast to 
the much faster, multi-peak variability observed during the prompt 
GRB. 

This Letter proposes an alternative possibility that the flares 
are produced by the long-lived reverse shock when it crosses the 
tail of the GRB ejecta. Our model does not require a long-lived cen¬ 
tral engine; it requires that the Lorentz factor in the end of the GRB 
explosion is significantly reduced, from L ~ 100 to L ~ 10, which 
leads to the formation o f an extended tail te.g. lUhm & BeloborodovI 
l2007l : l^net et al.l2007h . As described in Section|2] dense shells are 
expected to form in the tail of the expanding ejecta before the pas¬ 
sage of the reverse shock. In Section[^we estimate the bolometric 
light-curve of emission produced by the encounter of the reverse 
shock with the dense shell, and find that it resembles the observed 
X-ray flares. Comparison of the model with observations and future 
possible ways of its development are discussed in Section|3 


2 DENSE SHELLS AFTER INTERNAL SHOCKS 


The observed variability of GRB emission 

Fishman & MeesanI 

1 19951; iBeloborodov et al.ll2000l;lGuidorzi et al. 

|2012|) suggests that 


the Lorentz factor of the relativistic ejecta L fluctuates in a broad 
range of timescales, from a few milliseconds to several min¬ 
utes. In the expanding outflow, regions where the radial gradi¬ 
ent of Lorentz factor is negative (i.e. where L is decreasing out¬ 
wards) are progressively compressed. If the ejecta is not magnet- 
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ically dominated, internal shocks eventually form and propagate. 
The internal shocks can impact the prompt y-ray emission mainly 
in two ways: (i) in the subphotospheric (optically thick) region, 
shock heating can offset the adiabatic cooling of radiation and 
change its spectrum from thermal to the observed Band-type shape 
through t he Comptonization process and additional synchrotron 
emission l Meszaros & Reesll2000l : IPe’er et alJiTOOd : Iseloborod^ 
l20ld:IVurmetal.ll201lh : (ii) outside the photosphere, shocks can 


continue to produce synchrotron emissio n in the gamma-ray band 
if they efficiently accelerate particles i Rees & Meszarod 1 19941 : 


iKobavashi et '^ll997l:lDaigne & Mochkovitchlll998l)~ 

The outflow is likely initially inhomogeneous, and internal 
shocks can amplify the density variations, especially if the shocks 
are radiatively efficient. More importantly, the outflow tends to de¬ 
velop a few massive dense shells long after the shocks. At time to 
when the internal shock phase ends, the fluctuations of Lorentz fac¬ 
tor r have been damped and T acquires a monotonic radial profile. 
Its slope dr/dr > 0 is significantly inhomogeneous, as illustrated 
below by a hydrodynamical simulation. The subsequent ballistic 
expansion of the outflow at f » to generates a large contrast in 
density, as can be seen from the following estimate. 

Consider a shell of width Ash with the Lorentz factor variation 
across the shell dT = {dTIdr) Ash. As the shell expands ballistically, 
its width evolves as 


Ash(o ^ i + 

Ao r to 


( 1 ) 


where Ao is the shell width at to. The relevant Ao at the end of the in¬ 
ternal shock phase corresponds to the longest variability timescale 
observed in the burst, Aq ~ cfo, which is comparable to the burst 
duration, e.g. to ~ 1 - 10 s. In contrast, the time t given to ballis¬ 
tic expansion before the reverse shock crosses the outflow in our 
model is associated with the observed time of the X-ray flares, 
t~tf~ 10^ - 10“ s. 

The large ratio tjto implies that dT should strongly affect the 
density structure of the outflow. The parts of the outflow with dT ~ 
T will become much thicker than the parts with dT <K T, which will 
form dense shells. The corresponding density contrast is roughly 
given by 


pi 

-mm 

Pi 



( 2 ) 


An important feature of the post-internal-shock flows is the 
presence of a few “plateaus” in the profile of Tfr) or T(m) where 
m is the Lagrangian mass coordinate in the flow. These plateaus 
contain a large mass and eventually form dense massive shells be¬ 
cause of their small dT. We have observed and studied this ef¬ 
fect in outflows with various initial r(m) using detailed hydro- 
dynamical simulations. We found that that the plateaus typically 
have a small dT/T of a few per cent, even in the case of adia¬ 
batic shocks (an example is shown in Figure [T}. The simulations 
have been performed using a one-dimensional (spherically sym¬ 
metric) hydrodynamic code with a second-or der HLLC approxi- 


mate Riemann solver i 

Mienone & Bodoll2005b on a movina mesh 

dDuffell & MacFadven 

2011). 


Similar (but exactly flat) plateaus are produced by sim¬ 
plified simulations in which the outflow is represented by a 
large number of discrete shel ls that interact by direct collisions 
JPaigne & Mochkovitchll99^ . In this model, the ejecta is assumed 
to be cold and all pressure waves are neglected. This simplification 
is not crucial as it gives results similar to the detailed hydrody¬ 
namical simulations (cf. Daigne & Mochkovitch 2000). Below we 



Figure 1. Hydrodynamic simulation of a spherical outflow undergoing in¬ 
ternal shocks during its expansion. The flow is injected cold with a constant 
kinetic power during t„ = 10 s at a radius t?o = 0.1 ct„. The figure shows 
the Lorentz factor profile as a function of normalized Lagrangian mass co¬ 
ordinate, at three moments of time: the initial (injected) state (dashed line), 
during the internal shock phase at tiab = 800 s (dotted line), and after inter¬ 
nal shocks ended at fiaj, = 10^ s (solid line). At fiaj, = 10^ s, well after the 
internal shock phase, the outflow has cooled adiabatically and entered a bal¬ 
listic phase. As a result of internal shocks, two large portions of the ejecta 
have a uniform Lorentz factor with dT/T g 0.01. We have checked that 
adding a fast variability component in the outflow (leading to several gen¬ 
erations of internal shocks) does not significantly change the final Lorentz 
factor distribution. 


use the simplified model and extend the simulation to include the 
external blast wave and the reverse shock propagation through the 
outflow. 


3 FLARES: A SIMPLIFIED MODEL 

We used the simplified hydrodynamical simulations to follow the 
reverse shock (RS) propagation through the outflow and its inter¬ 
action with the dense shells. A sample model is shown in Figurej^ 
In this example, the explosion is assumed to produce an outflow 
of variable Lorentz factor distribution and duration /„ = 10 s. The 
injected kinetic power £k = 10^'^ erg s“* is assumed to be constant. 
Figure 2 shows how the initial distribution of F is changed after the 
internal shocks. Most of the outflow mass concentrates into three 
shells with Lorentz factors F x 240, 120 and 30 carrying about 10, 
20, and 40 per cent of the total mass. The last 30 per cent of the 
outflow have not been atfected by internal shocks. 

We followed the outflow to much larger radii where the RS 
eventually crosses all three shells. The dynamics of the reverse 
shock depends on the external density. We have calculated two 
cases: a uniform external medium p = const and a wind-like 
medium p = Ar^^, where A is a constant which depends on the 
mass-loss rate of the GRB progenitor. Its ty pical range for g alac- 
tic Wolf-Rayet stars is 10‘* - 10*^ g cm“' icrowthen 120071) . Re¬ 
cent modelling of optical and GeV flashes in GRB s suggest that 
A ~ 10*' g cm * is typical for GRB progenitors dHascoet et al.l 
l2014al : Hascoet, Vurm & Beloborodov, to be submitted). 
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prompt afterglow 




Figure 2. Left: Lorentz factor distribution as a function of Lagrangian mass coordinate; the front edge of the outflow is at m = 0. The initial distribution r(m) 
is shown hy the dashed curve, and the final distribution (after internal shocks) is shown by the solid curve. Right: bolometric light-curves (assuming a radiative 
efficiency ee = 0.1) calculated for different ambient media (dotted lines). The wind density p is parametrized hyAn as follows: pR“ = 10'*An gcm^'.Inthe 
uniform medium case, n is the density in cm“^. From top to bottom: wind medium with An = 1, uniform medium with n = 1 (Lboi has been divided by 10^), 
and the case of a “naked” burst (no ambient medium; Lbol has been divided by 10^). The solid lines show the same light-curves for the case of an anisotropic 
emission, which is beamed within an angle 6beam ~ 60° in the plasma frame. In the case of the naked burst, the three eaidy pulses (coiTesponding to the three 
pairs of internal shocks that initially propagate within the ejecta) are followed by high-latitude (off-axis) emission. 


3.1 Bolometric light curve 

The simulation gives the power dissipated in the RS and its bolo¬ 
metric luminosity with two assumptions: (i) that a fraction Cc ~ 0.1 
of the dissipated energy is injected in shock-accelerated electrons, 
and (ii) that these electrons promptly radiate their energy (“fast¬ 
cooling” regime). To find the bolometric light curve received by 
a distant observer we transform the emission from the flow frame 
to the observer frame and take into account the spherical curva¬ 
ture of the emitting shells. The radiation received by the observer 
from the RS is shown in Figure 2. In this figure, we also included 
the prompt radiation produced by the internal shocks, assuming the 
same Ce = 0 . 1 . 

When the reverse shock goes through the dense shells there 
is a sudden rise of dissipated power producing a flare in the light 
curve. The immediate rise of luminosity to the peak of the flare in 
Figure 2 is an artefact of the simplified model, which does not re¬ 
solve the propagation of the reverse shock through the dense shell. 
An accurate hydrodynamical model would give a slower but still a 
steep rise; it is discussed and estimated in Section|4]below. 

In our example with three dense shells, two flares are pro¬ 
duced; the crossing of the first shell (at F 240) corresponds to 
the initial rise of dissipated power from the reverse shock. Using 
the approximate expression for the flare radius Rf ~ F^ctf, one can 
estimate the observed time of the flare, 

= ^84o( —F-'/'s (uniform), (3) 

\ An n nip J \ n J 

= . ^3 r 4 ~ (t^) ^ (wind). (4) 

AnAc^T"^ \^ii/ 

Here F is the Lorentz factor of the massive shell (F 2 = F/100), E is 


the kinetic energy of the ejecta that has already crossed the reverse 
shock (£54 = F/IO^"* erg), n is the density of the uniform medium 
in cm“^ and An is the density parameter of the wind medium in 
units of 10 " g cm“'. One can see that the flare can be significantly 
delayed and this delay is sensitive to the Lorentz factor of the mas¬ 
sive shell. 

We also observe that the flares produced by this mechanism 
share two key features with the observed X-ray flares: (i) after the 
flare the light curve returns to the pre-flare decaying afterglow and 
(ii) the flare duration Aff is proportional to its time of occurrence tf. 

The simulation shown in Figure 2 assumes fast cooling at all 
times, which gives a constant radiative efficiency equal to e^- The 
flux enhancement during the flare can be even higher if before and 
after the flare the emission is in the slow-cooling regime. Then the 
radiative efficiency can be significantly increased during the flare. 


3.2 Effect of anisotropic emission 

In the simplest model, which assumes isotropic emission in the 
fluid frame, the characteristic temporal width of the flare is Atf ~ r 
where 


Then the ratio Atf/tf is close to unity, greater than the typical ob¬ 
served ratio Atf/tf ~ 0.1 - 0.3. 

Significant anisotropy is, however, expected and supported 
by the fast luminosity v ariations observed in GRB afterglows 
( [Beloborodov et all 1201 ih . The effect of anisotropy on the flare 
light curve is demonstrated by the following simple model. 
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For the moment, let us picture the encounter of the reverse 
shock with the dense shell as an instantaneous flash of energy Ef at 
a radius R. The bolometric light cu rve received from such a flash is 
L by jBeloborodov et alfcoilh . 


given I 


L(f) = 


2EfA{e) 




(t > tf). 


( 6 ) 


where 


cos 6 = 


T-(t- tf) 
T + it-lfY 


(7) 


6 is the angle with respect to the radial direction in the fluid 
frame, and A(6) describes the angular distribution of emission; for 
isotropic emission A(d) = 1. The magnetic field in GRB shocks is 
in the shock plane, and their synchrotron radiation is anisotropic 
even when the emitting electrons are isotropic; in this case A(9) = 
(3/4)(l +cos^ 9). Anisotropy of the electron distribution can further 
enhance the anisotropy of the emitted radiation, which significantly 
shorte ns the duration of the flare (see Figure 3 in iBeloborodov et aP 
l201lh . leading to a characteristic duration Atf <K tf. The shorter du¬ 
ration corresponds to a steeper decay of luminosity after the peak, 
with initial temporal index a ~ 3tf/Atf. 

As an illustration, in Fig. we show the model where the 
emission is moderately limb-darkened in the plasma frame, so 
that it is beamed within an angle 0beam ~ 60° (see Equation 9 in 
iBeloborodov et al1l201 ih . Then Atf ~ 0.3ff and the resulting light 
curves resemble the observed X-ray flares. 


4 SHOCK REFLECTION FROM THE DENSE SHELL 


In the simplified model of Section 3, the RS crossing of the thin 
dense shell instantaneously dissipates part of the shell kinetic en¬ 
ergy, and the observed light curve of the flash is determined by the 
curvature effect. In this section, we describe a more realistic pic¬ 
ture of the RS interaction with the dense shell. We also estimate the 
timescales for the rise and duration of the dissipation event. 

Two main effects should be taken into account: (i) The RS is 
greatly slowed down as it enters the dense gas, and it takes a finite 
time (estimated below) for the RS to cross the shell, (ii) The RS by 
itself does not efficiently decelerate the shell. The shell continues 
to plow through the ejecta that was previously accumulated in the 
blast wave. This ploughing launches a new internal shock inside 
the blast wave (Figure O. Effectively, the RS splits in two when it 
enters the dense shell. One (slow) shock continues to propagate into 
the dense shell and the other (fast) shock is reflected and launched 
back into the blast wave material. The fast shock dissipates most of 
the kinetic energy of the massive shell measured in the rest-frame 
of the blast wave. Below we call it ‘internal forward shock’ (lES). 
It is the lES that emits the X-ray flare. 

The velocity of RS propagation through the dense shell, Yrs = 
vrs/c ^ L may be estimated using the approximate pressure bal¬ 
ance between the RS and IFS. Let p 2 be the proper density of the 
shell just before it is reached by the RS. Assuming that the shell is 
cold and weakly magnetized, its pressure behind the RS jumps to 



(8) 


Now, let us estimate the pressure behind the lES, pips- The dense 
shell acts like a wall and the Lorentz factor downstream of the lES 
is almost equal to the shell Lorentz factor Fej. The relative Lorentz 
factor of the upstream and downstream is 



blast wave of thickness ^ 


Figure 3. Schematic illustration of the reverse shock (RS) encounter with 
a thin massive shell. The RS becomes slow as it enters the dense material 
and a new ’internal forward shock’ (IFS) is created, which propagates away 
from the dense shell into the blast wave material previously shocked by the 
RS. The figure also shows the forwai'd shock (FS) in the external medium 
and the contact discontinuity (CD) that separates the ejecta from the exter¬ 
nal medium inside the blast wave. 




(9) 


where Fb„ is the Lorentz factor of the blast-wave material upstream 
of the IFS, and we have used F^j » I and F^w » 1. Gas pres¬ 
sure m ay be written as p » (pi</3)(l -l- c^/n) jBeloborodov & Uhml 
l2006ll . where p is the density and u is the energy per unit rest mass. 
The gas in the downstream of IFS has been shocked twice: first 
by the RS and then by the IFS, both times with the same relative 
Lorentz factor Lrel- Therefore, u x F^^jC^ and p (TF^ij^pi, where 
Pi is the outflow proper density prior to the crossing by the RS and 
IFS. In the estimate for p we used the compression factor of 4Frei 
for a strong adiabatic shock. This gives 

Fifs ^ y (rtei - (10) 

The approximate pressure balance prs ~ Pifs gives an estimate for 
/Irs. 


8 

As- 3 


P2 


1/2 


( 11 ) 


Using the obtained j8rs, we estimate the time it takes the RS 
to cross the dense shell, /cross- The shell thickness is comparable to 


A3h = t/|^jAb„, (12) 

where Abw ^/r^w ~ 1® ll*® blast wave thickness, and rj = 

msb/mbw :£ 1 is the ratio between the mass of the incoming dense 
shell, and the outflow mass already accumulated in the blast 
wave, mbw. Then we find 
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Typical expected pilpi > 30 implies a short fcross- 

The shell Lorentz factor is reduced from Tej to T^w when it 
sweeps up the rest-mass 


'"sh 



(14) 


where we took into account that the gas ahead of the shell has been 
heated by the RS, and so its inertial mass is increased by the factor 
of Tpei. Assuming a relativistic IPS speed ~ c, this gives the follow¬ 
ing estimate for the observed duration of the shell deceleration: 


Abw /Tbw) ^ /T^w 

iflhvj r y Tej j \ Pej 


(15) 


The IPS may be mildly relativistic, Prei ~ 2, if the flares occur 
long after the formation of the ballistic outflow (t » tg). If the blast 
wave dynamics between the flares is approximated as self-similar, 
one finds Pcj/Pbw = V2 and 2 in the case of the wind and uniform 
external medium respectively. 

The time tcross is associated with the rise of the energy injec¬ 
tion event and f^iss determines its duration. The above results give 
typical fcross fdiss t. Por short ?diss the observed duration of the 
flare will be determined by the curvature (and anisotropy) effects, 
as described in Section 3.2. Here we note an additional effect which 
reduces the observed duration of the flare. While the pre-flare after¬ 
glow is emitted by gas moving with Pi,„, the flare is emitted by gas 
moving with Pej > Pb* (the gas behind the IPS). Therefore, even 
without the anisotropy effect, the curvature timescale for the flare 
is compressed by the factor of Atf/lf ~ (Pbw/Pej)^- 


5 CONCLUSIONS 

This Letter suggests that GRB outflows naturally develop a few 
dense shells with relatively low Lorentz factors, which leads to a 
new explanation for the X-ray flares. The dense shells give sudden 
energy injections to the GRB afterglow. We examined the dynamics 
of the injection event and showed that it generates a new internal 
shock behind the blast wave, which produces a short-lived flare. 
This mechanism provides a natural explanation for the fast rise of 
the X-ray flares and their observed duration Istfltf ~ 0.1 - 0.3. 

The proposed scenario assumes two radiative properties of 
the shock: (i) a fraction of the shock energy is transferred to non- 
thermal accelerated electrons capable of emitting synchrotron X- 
rays, and (ii) the electrons radiate X-rays in the fast cooling regime. 
Detailed modeling of hydrodynamics and radiative properties of the 
flares are deferred to a future work. As described in Section 4, the 
flare is emitted by gas that has been shocked twice and has high 
density and pressure. This gas can also possess strong magnetic 
fields and can radiate more efficiently compared with the pre-flare 
afterglow. This will help to satisfy the conditions (i) and (ii). 
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